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Abstract 

We discuss consequences of finite plasma formation time at RHIC. Nuclear modification factor, 
azimuthal assymetry, di-jet correlations can quantitatively be described by particle production in 
the early stage of the nuclear collision. A possible impact on interpretation of the non-photonic 
electrons and data is also considered. 
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The existence of a long Quark Gluon Plasma (QGP) formation time Q], T=2-3 fm/c, 
may have a strong inpact on some physical observables in the experiment. Direct physical 
interpretation of these observables could be missleading. If jets or other hard scattering 
processes do not suffer from strong energy loss or absorption during formation time, 
there should be a significant corona effect. All processes in the corona region will not be 
influenced by the produced dense medium. The relative contribution of corona decreases 
with centrality, but could account for at least 20% of hard processes in the most central 
collisions. It means, for example, that at each centrality class, the nuclear modification 
factor Raa for a particle produced in a hard scattering proccess can't be smaller than the 
corona contribution. Of caurse, this is valid if an additional absorption, like normal nuclear 
absorption, does not takes place during these first 2-3 fm/c. Many physical observables will 
be influenced by contribution from the corona region. A list of possible cosequences follows. 
We used a simple Monte Carlo simulation of nucleus-nucleus collisions based on Glauber 
approach with a Woods-Saxon nuclear density distribution jlj|. 

1. The PHENIX collaboration at Relativistic Heavy Ion Collider, RHIC, showed a pre- 
liminary result for the Raa dependence on azimiithal angle relative to the reaction plane in 

Fig. [T] As expected, Raa is larger in 



Au+Au collisions at v /s/v r /v=200 GeV at RHIC 
the reaction plane and smaller out of the plane. The most interesting feature of this result 
is that in the event centrality class 50-60%, the in-plane Raa becomes close to one. That 
is, no absorption at all is seen for high p t pions. At the same time a significant particle 
absorption is seen in the out-of-plane Raa- Here Raa is a nuclear modification factor which 
is defined as 

, x _ d 2 N A+A /dp T drj 
AA[PT) ((Nu^/a^dW^/dpTdr)' 1 j 

where (Nb inary ) is a number of binary nucleon-nucleon collisions at a particular centrality 
class. In the 50-60% event centrality class, the amount of nuclear matter is still significant in 
all directions. It's puzzling that a high momentum pion can "punch through" the interaction 
zone in the in-plane direction but may be stopped when traveling in the other direction. In 
the paper we attempted to explain numerically this feature by a QGP formation time 
of T = 2.3 fm/c. During this time, for certrality class 50-60%, in-plane produced jets can 
freely leave the interaction zone. From other side, jets produced out-plane, where the size of 
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interaction region is larger, have not enough time to escape and will be partially absorbed 
by the formed dense medium. 

2. Obvious consequences of such a scenario are: a) all pions above 5-6 GeV/ c are produced 
by parton fragmentation from the corona region, which is formed between 2 to 3 fm/c after 
the beginning of the collision; b) combining this with the experimental fact that the hadron 
to pion ratio, enhanced at low momenta, returns back to its vacuum value at pt above 5 
GeV/c we can conclude that all high p t hadrons are produced from corona region; c) 
Raa versus p t for pions and hadrons at high transverse momenta will be flat ^J. 

3. The model parameter T = 2.3 fm/c was adjusted to describe a small subset of Au+Au 
data at centrality 50-60%. However, with this value of T we can nicely explain the Raa 
behaviour not only for all centrality classes of Au+Au collisions, but for Cu+Cu interactions 
too, Fig. U 

4. We also estimate another experimental variable, v%. It describes the elliptic shape 
of the nucleus-nucleus collision in azimuthal angle as dN/d0=N(l+2u2Cos(2(/>)). The long 
formation time naturally brings the result of non-zero vi for high p t pions Fig. El The 
value of t>2 reaches 11-12% in mid-central collisions of gold nuclei with Woods-Saxon density 
distributions. At each centrality class i>2 will be constant at high transverse momenta. 

5. Two-particle angular correlations of hadrons close to the high p t "trigger" particle 
or opposite to the "trigger" particle directions allow to investigate of the near-side and the 
away-side jets properties, respectively. Following the argumentation that all particles at high 
p t are produced from corona region, we can explain the weak dependence of properties of 
the near-side jets j]| with centrality in Au+Au collisions. In the cited paper, some increase 
was observed in the near-side particle yield only for low momentum associated particles. 

6. STAR collaboration published the results on Iaa in Au+Au collisions - the ratio of 
away-side yield per trigger particle to the similar value from p+p collisions ^]. We calculate 
Iaa within our model with T = 2.3 fm/c for different centralities. The increase of corona 
thickness by 2.3 fm keeps the value of Iaa at the level close to Raa, Fig. EJ We want to 
emphasize that these are two different variables, Iaa is related to the yield per trigger 
particle. The trigger particle should "survive" from absorption in a first place. In this 
calculation we assume that away-side jet distribution has a Gaussian shape with the width 
parameter of 0.35 radians. 

7. At Quark Matter 2005 conference, August 4-9, in Budapest, Hungary, STAR col- 
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laboration claims that away-side jet re — appears at high trigger and associated particles 
momenta While changing the momentum range of the associated particle, the jet width 
will change too - for higher momentum particles the jet width should be smaller. We check 
the STAR result by changing the width of simulated away-side jet, Fig. For all values of 
width, Iaa is still significant even for most central events. We argue that in previous analysis 
at lower momentum range, the away-side jet was wider and it was not easy to reconstract 
its remnants on the level of 20-25%. Besides this, at low associated particle momentum, the 
responce of the high density medium on the absorbed away-side jet changes dramatically in 
shape and in yield |lfl| . 

8. Discussions in the previous parts 2, 6 and 7 lead to the conclusion that, probably, 
even at particle momenta 5-10 GeV/c, we don't see "punch-through" jets, these are di-jets 
from corona region. This is why width of the away-side jet does not change with centrality 
and is the same as in p+p collisions [9j. We investigate the away-side jet shape in our 
model. Di-jet forms kind of a tangential emission from the surface of interaction region. If 
in experiment the trigger particle radiates at particular direction, di-jets can survive along 
the two tangents, from the edges of the collision zone. If the collision zone is small and 
away-jet is wide we may see jet splitting too. In Fig. [Hjwe plot the results of the simulation 
of di-jets in mid-central 30-35% Au+Au events. The away-side jet originally had a Gaussian 
shape with the width of 0.75 radians. The dashed histograms represent jets survived at the 
two tangential directions. We see slight shifts from the original direction on the order of 
25% of the original width. The solid line histogram represents the sum of these two subsets 
and its fit by the Gaussian function. The width parameter sigma of this fit, 0.77 radians, 
is very close to the original width. The observed features do not change significantly with 
centrality. We can conclude that, indeed, the shape and the width of away-side jet are not 
changed much. There is no visible back jet splitting produced in the corona region, thus 
it can not be an explanation of the two peak structures seen in the experiment [lOj at low 
momentum of the associate particle. 

9. In Fig. Owe present PHENIX preliminary results on non-photonic electrons at high 



momenta shown at Quark Matter 2005 11 1 . We also plot the lower limit of Raa determined 
by particle production from the corona region during the formation time, solid line. Exper- 
imental points sit exactly on the line. These electrons originated from D-mesons, charmed 
baryons and, probably, from bottom quark decays. The conclusion of Fig. would be that 
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at high p t the charm quarks are completely absorbed by dense medium, similar to the light 
quarks. 

10. In the case of a strong charm quark absorption, following the argument in our previous 
part 4, one can expect a significant v 2 for non-photonic electrons. Experiments indeed found 
non-zero v 2 jlllQl- Qualitativety, it is possble to descibe also the momentum dependence 
of t>2. Leptonic modes of D- meson decay are mostly with 3- and 4-particles in the final 
state. The low momentum electrons should be produced isotropically from the low energy 
D-mesons. With increasing energy of the meson, decay products will start to follow the 
direction of deca ying meson, so electron v-i will increase too (this argument was also stressed 
in the reference |l2j), and finally will saturate for high p t at the values presented in Fig. 01 
It would be mistaken to say naively "charm flows". This statment needs clearefication and 
an additional experimental investigation. 

11. We expect many interesting features from J/ip production in nuceus- nucleus collisions. 
J/if> is a closed system of two heavy charm quarks. In particular interest is the interation 
of such a system with the surrounding medium. Previously, a very strong J/ip suppression 
was observed at lower energy nucleus-nucleus collision s |l4l Il5l | . New PHENIX preliminary 
data show the suppression at about similar strength [13j , Fig. |H| In our consideration we 
devide nucleus-nucleus collision in to two stages: before plasma formation and after. As it 
was already demonstrated, hard scattering processes, like jet production, do not suffer from 
energy loss or absorption during plasma formation time. The matter is transparent for them 
at this stage. In contrast, during the actual "plasma time", the medium looks completely 
"black" for jets. 

Part of J/ips could be absorbed during early stage of the collision. In this case the medium 
will be "grey" for them. At the same time, other models predict incomplete suppression of 
3 /if) in QGP phase, which means also that plasma is "grey" . We try to demonstrate how the 
distinction of plasma formation time, QGP itself and different suppression processes lead to 
a few major scenario, Fig. |H1 Scenario (a), white and black: no absorption at the early stage, 
total absorption in the QGP - solid line in Fig. |H1 It this case everything is similar to the jet 
absorption. Scenario (b), grey and black: during the plasma formation time J /if) has normal 
nuclear absorption but suffers from total absorption in the QGP, lower dashed line. This 
line goes obviously below the solid curve and was obtained by applying the results of R.Vogt 
calculations |lf| to the corona region only. In other words, it is the solid line result from 
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the ) times the original curve from R.Vogt. This would be an overestimation of the 

normal nuclear absorption because we did not take into account that the initial stage lasts 
2.3 fm/c only. Scenario (c), white and grey: top dotted line, where we assume that there 
is no absorption during the initial stage and some suppression in the QGP stage. Here we 
take just one model with QGP absorption by Kostyuk et al. [l^ and applied to the second, 
QGP, phase only. 

It is worthwhile to mention that the original calculations jltl Il7| for J/^ we used in cases 
(b) and (c) significantly overestimate and underestimate Raa, respectively jl3{ ]. 
By listing such naive speculations we want to stress the importance of the different stages 
of the nuclear collision, including plasma formation time, for J /if; production. 

In conclusion, we want to emphasize that we propose not "yet another" simple model with 
just one parameter, which actually, may increase a little bit if the longitudinal expantion 
will be taken into account. We do not pretend to describe in many details all features of 
the nucleus-nucleus collision at RHIC. We present an orthogonal point of view on some 
experimental data. This view could be useful for better understanding of properties of a 
new matter produced at RHIC. Perhapes, it will bring clear vision or some new ideas on 
formation process and properties of a quark-gluon phase. 

We tried to explain some of the features seen in the experiment, but also bring two 
major questions: where this formation time comes from, why models with fractional parton 
energy loss do not work well or why energy loss is very strong? The answers may come 
soon: interesting new approaches such as the complex plasma model or the radically 
new picture of QCD with polymer chains have been proposed recently. 
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For all experimental data with the stamp "preliminary" we scaned avalible electronic version 
of the plots by program xyscan developed by Thomas S. Ullrich. 
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FIG. 1: PHENIX preliminary Raa results for tt° at momenta 5-6 GeV/c versus angle <j) relative to 
the reaction plane for different centralities in Au+Au collisions at RHIC 3]. Points are experimental 
data, thin lines show systematic errors from the reaction plain resolution and vertical bars in the 
middle show the reaction plane avaraged Raa value and its error. Black horizontal bars are 
predictions of our model. 
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FIG. 2: Calculated Raa for Au+Au, dashed curve, and Cu+Cu, solid curve, collisons versus the 
number of participant nucleons, N par t. The circles are experimental data for 7r° yield for Au+Au 
intergrated for pt><i GeV/c Q|. The triangles are data for Cu+Cu intergrated for pt>7 GeV/c 
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FIG. 3: Calculated ellipticity parameter v% for Au+Au collisons, solid line, versus the number of 
participant nucleons, N part . Data for 7T° with error bars are: squares for 5-7 GeV/c, cicles for 4.59 
GeV/c. PHENIX preliminary data HQ. 
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FIG. 4: Calculated Iaa for Au+Au, the curve, versus the number of participant nucleons,iVp art . 
The width (sigma) of away-side jet is 0.35 radians. The cicles are experimental data from 
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FIG. 5: Calculated Iaa for Au+Au versus the number of participant nucleons, N part , for different 
width of away-side jet - the numbers in radians next to the lines. Results for Cu+Cu at 200 GeV 
and Au+Au collisions at 62.4 GeV are also shown. 
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FIG. 6: Angular distributions within the away-side jet in the case of jet production from the 
corona region with strong jet absorption in QGP core. The original away-side jet had a Gaussian 
shape with the width parameter of 0.75 radians. Dashed histograms represent the jet distribution 
from two different sides of the collision zone. The solid line histogram is the sum of these two 
subsets. Numbers in the box represent the result of the solid line histogram fit by the Gaussian 
function. 
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FIG. 7: Raa for non-photonic electrons in Au+Au versus the number of participant nucleons, 



Q: 



N par t. The PHENIX preliminary experimental data are from circles and squares are for 

Pi=4.25 GeV/c and 4.75 GeV/c, respectively. Shown are statistical errors only. Solid line - the 
result of our calculation which puts the geometrical limit in case of strong charm quark absorption. 
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FIG. 8: Raa for J/i/j in Au+Au versus the number of participant nucleons, N par t. The PHENIX 
preliminary experimental data are from cicles - at midrapidity, squares - at rapidity y = 

1.2 — 2.2. Shown are statistical errors only. For explanations of the different lines, see discussion 
in the text. 
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